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ABSTRACT 

Thermonuclear explosions may arise in binary star systems in which a carbon-oxygen (CO) 
white dwarf (WD) accretes helium-rich material from a companion star. If the accretion rate 
allows a sufficiently large mass of helium to accumulate prior to ignition of nuclear burning, 
the helium surface layer may detonate, giving rise to an astrophysical transient. Detonation of 
the accreted helium layer generates shock waves that propagate into the underlying CO WD. 
This might directly ignite a detonation of the CO WD at its surface (an edge-lit secondary 
detonation) or compress the core of the WD sufficiently to trigger a CO detonation near the 
centre. If either of these ignition mechanisms works, the two detonations (helium and CO) 
can then release sufficient energy to completely unbind the WD. These "double-detonation" 
scenarios for thermonuclear explosion of WDs have previously been investigated as a potential 
channel for the production of type la supernovae from WDs of around one solar mass. Here 
we extend our 2D studies of the double-detonation model to significantly less massive CO 
WDs, the explosion of which could produce fainter, more rapidly evolving transients. We 
investigate the feasibility of triggering a secondary core detonation by shock convergence 
in low-mass CO WDs and the observable consequences of such a detonation. Our results 
suggest that core detonation is probable, even for the lowest CO core masses that are likely to 
be realized in nature. To quantify the observable signatures of core detonation, we compute 
spectra and light curves for models in which either an edge-lit or compression-triggered CO 
detonation is assumed to occur. We compare these to synthetic observables for models in 
which no CO detonation was allowed to occur. If significant shock compression of the CO 
WD occurs prior to detonation, explosion of the CO WD can produce a sufficiently large 
mass of radioactive iron-group nuclei to significantly affect the light curves. In particular, this 
can lead to relatively slow post-maximum decline. If the secondary detonation is edge-lit, 
however, the CO WD explosion primarily yields intermediate-mass elements that affect the 
observables more subtly. In this case, near-infrared observations and detailed spectroscopic 
analysis would be needed to determine whether a core detonation occurred. We comment on 
the implications of our results for understanding peculiar astrophysical transients including 
SN 2002bj, SN 2010X and SN 2005E. 

Key words: hydrodynamics - radiative transfer - methods: numerical - binaries: close - 
supernovae: general - white dwarfs 



1 INTRODUCTION 

Type la supernovae (SNe la) are understood to result from the ther- 
monucle ar disruption of a carbon-oxyge n (CO) white dwarf (WD) 
star (e.g. lHillebrandt & Niem ever 2000). One possible mechanism 
for igniting such an explosion can occur in binary systems in 
which a primary CO WD accretes He-rich material from a donor 
star. When a sufficiently large surface He layer is accreted, it is 
expected to ignite expl osively leading to detonation of the ac- 
creted He layer (see e.g.lNomotdll980 L il982l : IWooslev et al J 198(1 
IWooslev & Weave J 19941) . Detonation of the He layer can then lead 



to a secondary detonation of the core, either by directly igniting 
the CO fuel near to the interface with the overlying He layer (see 
e.g. lNomotoll982l ; lLivne & Glasneilll990h or due to compre ssional 
heatin g of the core by inward propagating shocks (see e.g. iLivnd 
Il990h . The consequence of this "double-detonation" model is the 
incineration of the CO WD and its He-rich outer layer, leading to 
an explosion in which the primary star is completely destroyed. 

However, the question of whether the secondary detona- 
tion forms is challenging owing to the wide range of rele- 
vant length-scales that must be resolved if it is to be sim- 
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ulated dSeitenzahl etal] l2009al lbb. To date only a few multi- 
dimensional studies have been made of the double-detonation sce- 

r 



nario for a handful of progenitor mod els ( Livne & Glasne iill99C 



Dgani & Livio 1990; Livne & Gl asnedl991l;lLivne & Arne ttl 19951; 



Benzj 
2001 ; 



19971; iGarcia-Senz et alj|l999t iForcada et alj|200d : lForcad3 
Finketalj|2007l . l2010h . If the secondary core detonation 
does not occur, the result of He ignition is quite diff erent from the 
double -deto nation SN la model - as described by iBildsten et al.l 
d2007l) and IShen & Bildstenl d2009h . explosive burning of an ac- 
creted He-layer alone will lead to a thermonuclear transient that 
is roughly ten times fainter and evolves significantly faster than 
a SN la. Dubbed "point-la" (hereafter p-Ia), this class of explo- 
sion is readily accessible to observation by the current genera- 
tion of transient surveys (e.g. the Panoramic Survey Telescope & 
Rapid Response System [Pan- STARRS J], the Palomar Transient 
Factory [PTF0 and planned wide- field surveys by instruments such 
as SkyMapper [Keller et al. 2007] and t he Lar ge Synoptic Survey 
Telescope ILSST Science Collaboration! l2009ll ). Indeed, transient 
events with some similarities to the predicted properties of p-Ia 
explosions have already bee n reported and modelled in the con- 
text of the p-Ia scenario (e.g.lFolev et alj|2 009; Pere ts et al ] |20ld ; 
Poznanski et"ai]|2010l ; fkasliwal et alj|20ld ; IWaldman et alj|201 it 
Sullivan etalj|20 111) . It is to be expected that other similar events 



will be found and studied in the near future. 

In our previous studies dFink et allfeOlOtlKromer et alfcOlOh . 
we investigated the double-detonation scenario for systems with 
fairly massive CO cores (A/co > 0.81 Mg), Such cases are the 
most promising for yielding thermonuclear explosions as bright as 
normal SNe la. In those works, we focused on the possibility of 
core detonation triggered by converging shocks deep in the CO 
WD. We found that secondary co re detonation was very likely in 
all of the models we considered dFink etalj|2O10h . This implies 
that the p-Ia scenario should not be realised following He detona- 
tion in such systems unless some additional effect comes into play 
(see Section loTTV We now wish to extend our 2D studies of the 
double-detonation scenario to investigate systems with less mas- 
sive CO cores (e.g. Mco^0.6 Mq) and quantify the observable 
properties of double detonations for such systems. Compared to 
their more massive counterparts, these low-mass systems have two 
important differences. First, prior to any explosion, the central den- 
sity of the CO core will be lower. In principle, this might make 
it harder for the converging shocks to compress the centre suffi- 
ciently for a core detonation to occur. Previous 2D studies have 
already suggested that secondary detonatio ns can be produced fo r 
CO cores with masses as low as 0.55 Mq JLivne & Arnettll 19951) . 
Thus the first objective of our study is to extend our studies to even 
lower mass, around the minimum CO core mass that is expected to 
be realized in natu re (in the recent binary synthesis calculations of 
lRuiteretal]|201ll . the lowest CO core mass is ~0.45 M Q ). Sec- 
ond, even if a core detonation does occur, the low core density 
means it will produce little 56 Ni (see e.g. the lowest mass model of 
lLivne& Arnettll 19951 , which yields only 0.14 M Q of 56 Ni). In the 
limit of a very low 56 Ni-yield in the core, the radioactive products 
produced in the He-detonation may still play a dominant role in de- 
termining the explosion brightness and light curve evolution (as in 
the p-Ia scenario). Our second goal, therefore, will be to quantify 
the observable properties of explosions in which the core detonates 
but produces only a small mass ( < 0.15 Mq) of 56 Ni. 



We begin, in Section [2] by introducing the models adopted 
for this study and the suite of numerical simulations used to study 
them. In Section [3] we describe our simulations in which we in- 
vestigate shock convergence and detonation in low mass CO cores. 
We then present the results of alternative explosion simulations (in 
which either an edge-lit CO core detonation is assumed or in which 
it is assumed that no core detonation occurs; see Section^. In Sec- 
tion [5] we present the synthetic observables computed for all our 
simulations before discussing our results and drawing conclusions 
in Section[6] 



2 METHODS 

For this study we have performed sets of 2D numerical simula- 
tions that follow the explosion dynamics, nucleosynthesis and ra- 
diation transport of different explosions for two specific initial sys- 
tems. Here we describe the parameters of the initial systems and 
the means by which the numerical simulations were performed. 



2.1 Initial models 

As initial conditions for our explosion simulations, we adopt mod- 
els describing the state of the system immediately prior to He deto- 
nat ion. These are not based on evolutionary/accretion calculations 
(cf. lWooslev & Weaverlll994l ; IWooslev & Kasenl201lh but are ide- 
alized representations of low-mass CO WDs that have accreted a 
surface layer of He. Thus, in this work, it is an assumption that 
the system has evolved to reach conditions for He detonation - we 
only determine whether subsequent detonation of the CO core is 
probable and quantify its observable consequences. 

We have considered two sets of system parameters, which are 
described below. Population synthesis calculations suggest that 
He-rich accretion by CO WDs leading to helium shell detonation 
most commonly occurs i n binaries contain ing a CO WD primary 
and a He WD donor (e.g. lRuiter et alj201 llFI . However, that chan- 
nel is not predicted to yield significant numbers of systems with 
very low CO mass (^0.6 Mq). For the CO masses we consider 
a non-dege nerate He st ar companion is more probable (see figure 
3 of lRuiteretal.ll201ll) . The evolutionary path to a low-mass CO 
WD accreting from a He star involves one or more stable mass- 
transfer episodes and at least one common envelope phase leading 
to a tight binary (< 1 hr orbital period). Compared to binaries with 
He WD donors, systems that evolve to a final state involving a CO 
WD and a He star donor are rare r and the evolution ary timescale is 
usually relatively short (< 1 Gv; lRuiter et alj|201 ll) . Nevertheless, 
the He-star donor scenario is likely to be the most promising route 
to the low CO mass explosions that we will discuss. 

For our standard system (hereafter "S" model), we adopted a 
CO core mass of Alco = 0.58 Mq and a mass for the accreted 
He layer of MHe = 0.21 Mq, yielding a total mass of M tot = 
0.79 Mq. H ere, Mn e is close to th e mi nimum mass for detona tion 
suggested bv lBildsten et alj J2007l) and lShen & BildsterJ d2009h for 
our chosen value of Mco- This model naturally extends the study of 
iFink et al.1 < f2oToh into the regime of physically plausible low-mass 
systems that might be realised in systems where a primary CO WD 
accretes from a He-burning star (see Section [6j4) , Its masses are 



1 http://pan-starrs. ifa.hawaii.edu 

2 http://www.astro.caltech.edu/ptf 



3 Note, however, that the Ruiter et al. (2011) calculations assumed explo- 
sions occurred once a 0. 1 Mq layer of He builds up; the systems we discuss 
here require additional He accretion beyond that point. 
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very similar to the least ma s sive models considered in th e studies of 
IWooslev & Weaver! ( fl994h . iLivne & Arnettl lll995t) and lShen et al.l 
feoich . and also to model CQ.55HE.2 of lWaldman et al] ( 1201 lh . 

As a second case, we also considered an extremely low-mass 
system. This model (hereafter model "L") is designed to robustly 
bracket the low-mass end of the distribution of potential initial sys- 
tems. The adopted CO core mass (Mco = 0.45 Mq) lies at the 
lower boun dary of the distribut ion in the population synthesis cal- 
culations of iRuiter et alj d201ll) . We have also adopted a very low 
mass (Miie = 0.21 Mq) for the He layer when the explosion oc- 
curs. This is close to the most optimistic (i.e. lowest) estimate for 
th e mass of He needed f or detonation, following the arguments 
of IShen & Bildste n (2009fl Moreover, the evolutionary models of 



Table 1. Parameters defining the initial model. 



Parameter 



Model S Model L 



Wooslev & Kase nl (120111) imply that the conditions suggested by 
Shen & Bildstenl j2009h lead to He layer masses that are generally 



too small for detonation. Thus, our L-model likely lies outside the 
regime in which He detonation is probable. However, we include it 
as an important numerical experiment to test the limit of the double- 
detonation model - if He detonation in this system leads to CO 
core detonation, then it can be concluded that our method would 
predict secondary detonation for any combination of CO-core/He- 
layer mass for which He detonation is realist ic. Our L-model issim- 
ilar to the lowest mass model considered bv lWaldmanetalJfcOllh 
(C0.45HE.2). 

The initial models are set up in exac tly the same manner as 
the models described by iFink et all fcOlfjh by choosing appropri- 
ate values for the temperature and central density of the CO core 
(T c , p c ), and the temperature and density at the base of the He layer 
(21, Pb). These parameters are listed in TableQ] 

In our S-model we aimed at achieving a bright explosion 
and therefore assumed a cold shell (i.e. as dense as possible for 
a given mass). For degenerate matter, the exact value of T does 
not matter. Thus, we simply assumed a constant temperature T = 
5xl0 5 Kinthe whole WD. Following lshen & BildsteH j2009h and 
IWaldmanetailfeOllh . we adopted higher temperatures for our L- 
model (T = 1 x 10 7 K in the core and T = 2 x 10 s K at the 
base of the shell, decreasing adiabatically outwards). Higher tem- 
peratures reduce the density, leading to less complete burning and 
making it harder to trigger a secondary core detonation. Thus, this 
choice maintains the status of our L-model as a fairly extreme test 
for the plausibility of secondary detonation by shock compression. 

For simplicity, the core is assumed to consist of uniformly 
mixed 12 C and 16 (equal parts b y mass) and the s urface layer 
is assumed to be pure He. Following lFink et aT the density 

profile within the model is constructed by solving for hydrostatic 
equilibrium conditions with the adopted value of T c , p c , 21 and 
Pb us ing the same equation of state adopted by IFink et al] J2007I . 

[2oToh . 



2.2 Explosion simulations, nucleosynthesis and radiation 
transport 

In most respects, our s imulations were p erformed in the same way 
as those described by IFink et al] fcOld) . Therefore, we will only 
summarise the main points and highlight the small number of mod- 
ifications to the numerical implementation adopted. 



T C (K) 


5 x 10 5 


1 x 10 7 


p c (10 6 gcm- :i ) 


8.5 


3.81 


T„(K) 


5 x 10 5 


2 x 10 s 


p b (10 6 g cm" 3 ) 


1.3 


0.592 


M CO a (Mq) 


0.58 


0.45 


M He a (Mq) 


0.21 


0.21 


M tM a (Mq) 


0.79 


0.66 



From figure 5 oflsiien & Bildster] <2009l) . a minimum He-layer mass of 
~0.2 M© (for Mco = 0.45 Mq) is required, if one adopts the condition 
that the dynamical timescale is one tenth of the local heating timescale when 
dynamical burning sets in. 



a Note that the masses are not independent parameters but follow from T c , 
Pc, 7b an d Pb ar >d the assumption of hydrostatic equilibrium (see Fink et al. 
2010). 



All the simulations presented here were carried out in 2D with 
rotational symmetry about the z-axis. As in IFink et alj Jioicl) . 
begin by igniting a detonation at a single point in the He layer. As 
noted in Section [2~71 that such a He detonation ignites is a funda- 
mental assumption of our simulations. In all cases, we choose to 
ignite the He detonation at the base of the He layer on the positive 
z-axis. In the absence of evolutionary calculations prior to explo- 
sion, this is the simplest choice for He ignition but we note that it is 
disfavoura ble for edge-lit se condary detonations (see Section POt . 

As in lFink et al j teOlOl) . detonations in CO and He were mod- 
elled with a front tracking scheme using tabulated values for both 
detonation speeds (D) and energy release per unit mass (Q) be- 
hind the burning front (hereafter referred to as "detonation tables"). 
Since it takes place entirely in a low-density incomplete burning 
regime, the He detonation nucleosynthesis is very sensitive to the 
input parameters of the front-tracking scheme. The detonation ta- 
bles were therefore determined for each model of this work by 
applying the hydrodynami cs/post-pro c essing iteration scheme de- 
scribed in the appendix of IFink et alj | |2010|) . This time, however, 
the setup in the calibration runs was identical to that of the models, 
i.e., a detonation propagating laterally through the same He layers. 
To calibrate the He-detonation speed, the Rankine-Hugoniot jump 
conditions for detonations were solved for the minimum possible 
value of D, which corresponds to the flow velocity v^b of the final 
ash state being exactly sonic relative to the front. This procedure 
was repeated for every point on the tabulated grid of density val- 
ues. In the calculations we use the same equation of state as in the 
hydrodynamics code and the energy release from the previous it- 
eration step. This procedure leads to converged detonation speeds 
after around six iteration steps. The final detonation tables for our 
standard ("S") model are illustrated in FigureQ] It can be seen that 
complete He burning is never achieved for the densities present in 
the model and that this substantially reduces the detonation speed 
at low densities, compared t o the complete bur ning case. 

Unlike in the models of IFink et alj ( |2010|) . we suppressed any 
volume burning in the He layer before the arrival of the detona- 
tion wave. This provides a well-defined initial condition in the He 
layer (allowing the composition to change due to volume burning 
prior to detonation would not be self-consistent as we do not sim- 
ulate the evolution of the progenitor before the initiation of the He 
detonation). 

We performed detailed nucleosynthesis calc ulations for the 
explo sion models using a tracer particle method jTravaglio et alj 
120041) . Here, tracer particles are passively adverted in the hydro- 
dynamical simulation and used to record the thermodynamic tra- 
jectories of mass elements. These are then used as input to a 
post-processing step in which detailed isotopic yields are obtained 
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Figure 1. Detonation speed (Dtj; units of 10 9 cm s~ x ) relative to un- 
burnt matter (solid lines) and Q-values (energy release per unit mass in 
10 18 ergs g -1 ; dashed lines) for our S-model. The red curves indicate the 
final values obtained by our iterative calibration of the lateral detonation. 
The black curves show the values expected for a Chapman-Jouguet detona- 
tion if it were assumed that complete burning to 56 Ni occurs at all densities. 

from calculations with an extensive nucleosynthesis network (384 
species). We adopted an updated version of the REACLIB reaction 
rate library dRauscher & Thielemann! I2OOOL updated 2009) and a 
refined method with variable tracer masses dSeitenzahl et alj2010f) 
was applied. Variable tracer masses allow for a better spatial reso- 
lution at the edge of the CO core. 

The nucleosynthesis tracer particles are used to reconstruct the 
detailed abundances throughout the ejecta. This gives a complete, 
2D model for the structure of the ejecta at the final time of the 
hydrodynamical simulations (i.e. density and composition as func- 
tions of expansi on velocity i n the r- and z-directio ns). We then used 
the ARTIS code dSimll2007l: iKromer & Siml l2009) to compute syn- 
thetic light curves and spectra for the models. For all the radiative 
transfer simulations we used ou r set of atomic data extrac ted from 
CD23 of iKurucz & Belli d 19951) (see lKromer & Simll2009l) . but we 
expanded the range of ions included to I - VII for elements with 
atomic number 22 < Z < 28 to allow for higher ionization states 
that may be present at early times when the ejecta are hot. 



3 SECONDARY DETONATION BY CONVERGING 
SHOCKS; CSDD MODELS 

Our first simulat ions are the most similar to those described by 
iFink et al. I d201Ch - they study the shock convergence and poten- 
tial for formation of a secondary detonation via compression of the 
CO core in less-massive systems. We investigate this scenario for 
both our in itial systems (S an d L, see above). 

As in lFinketal.1 feoich . we simulate the propagation of the 
He detonation as it wraps around the CO core. The He detona- 
tion drives a shock front that propagates into the core leading to 
strong compression around a convergence point. Although the con- 
vergence point is off-centre in both models, it is less off-centre than 
in the models of IFink et al. I d2010h . This is a continuation of the 
tr end for ignition c loser to centre in less massive cores (see table 2 
of lFink et al.ll2010l) . 

Our first question is whether this compression leads to a large 
enough volume of sufficiently hot and dense material that a sec- 
ondary CO detonation could ignite. To assess this we compared the 



Table 2. Conditions at the hot spot in our converging shock double- 
detonation models. ij eI1 is the time at which critical conditions for ignition 
of the CO core detonation are reached while z- Kn / Rqo is the position of 
the hot spot (which lies on the 2>axis of the simulation) in units of the CO 
core radius. T Kn and p[ sn are the temperature and density at the hot spot at 
t = tjg„. A is the grid resolution. 



Parameter 


Model S 


Model L 


^ign (s) 


1.34 


1.81 


z ign (10 8 cm) 


-1.39 


-1.61 


Zign/RcO 


0.31 


0.31 


T ig „ (10 9 K) 


6.44 


4.68 


Pign (10 7 g cm" 3 ) 


18.0 


7.83 


A (10 6 cm) 


3.71 


4.94 



density and temperature reached in our sim ulation to critical tem- 
peratu res a nd densities for deto nation from iNiemever & Wooslevl 
dl997l) and iRopke et all d2007lfl From this comparison (see Ta- 
ble 13, we find that critical conditions for CO core detonation are 
robustly met in the simulation for both our initial systems - in fact, 
the peak tempera tures and densitie s are very similar to those of 
models 1 and 2 o f lFink et al.l d20ich . implying that core detonation 
is not significantly harder to achieve in the systems we consider 
here. Given the extreme properties of our model L, we therefore 
conclude that 2D converging shock simulations performed with our 
current approach will favour secondary core detonations for any 
physically plausible pair of CO core/He-layer mass. To study this 
further would require much higher resolution simulations that re- 
solve the critical volumes for detonation. 

Since our simulations of the shock convergence suggests that 
a detonation in the CO core is likely, we initiate a second deto- 
nation wave at the location of the hot spot in the CO core. This 
detonation sweeps over the whole CO core and releases sufficient 
energy to completely unbind the star. Hereafter, we will refer to 
the results of these simulations as our converging-shock double- 
detonation (CSDD) models (CSDD-S and CSDD-L for our two 
initial systems, respectively). 

The mass yields obtained from the nucleosynthesis post- 
processing of the CSDD models are tabulated in Table|3] The ejecta 
composition of model CSDD-S is illustrated in the top panels of 
Figure [2] which shows both the 2D distribution of mean atomic 
number and the detailed composition for a slice through the equa- 
torial plane of the model. 

For our CSDD-S model, a significant mass of radioactive nu- 
clei (specifically ~0.08 M of 56 Ni and 52 Fe) is produced by the 
detonation of the He layer. This yield of ra dioactive nuclei is simi - 
lar to that found bv lShen et all ( hoioh and IWaldman et alj d201ll) 
for models with comparable values of pb (spec ifically, our pat- 
tern of radioactive yields lies between those of Sh en et all 12010] 
for detonations of 0.2 and 0.3 Mq He-layers around 0.6 Mq CO 
cores). The decay of this material will power the early phases of 
the light curve and produce a transient that brightens rapidly, on 
a timescale of several days, as predicted by iBildsten et al.l d2007t) 
(see Section[5]l. In our model, however, the shock compression is 
sufficient to yield an even larger mass of 56 Ni (~0.1 M ) from 



5 Note that for the high densities reached in the shock convergence region, 
the critical volume for initiat ion of a detonati on is small compared to our 
grid resolution (see table 1 of [jink et al. 2007j). Therefore only temperature 
and density can be considered as detonation criteria in this case. 
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Figure 2. Composition of the explosion simulations for our S-model (CSDD-S, ELDD-S, HeD-S models, top to bottom). The left panels show the mean 
atomic number (Z) for each of the Lagrangian nucleosynthesis tracer particles at their final positions at the end of the hydrodynamical simulations (100 s 
after ignition of the He detonation). The models are symmetric under rotation about the 2-axis. Most of the dense CO core material in the HeD-S model (the 
encircled dark blue region around the origin ) remains bound and is not included in our radiative transfer simulations of the homologous ejecta. The right panels 
show the detailed composition for an equatorial slice through our simulations extrapolated to the homologous phase (the bound material in the HeD-S model 
is not included here). The black histograms show the total mass density (pt t) versus expansion velocity. The contributions from various important isotopes 
and elements are indicated by the coloured histograms in each plot; specifically, they show ptoi^;, where Xi is the mass fraction of the species in question. 



the core detonatio n. As in double-detonation models for more mas- 
sive CO WDs (e g.lWooslev & Weaver! 1994l ; lLivne & Arnettl 19951 ; 



iFink et al .l2007l ; IWooslev & Kasenl201lir this o6 Ni is concentrated 
at low velocities and surrounded by an envelope of intermedi- 
ate mass elements, predominantly silicon and sulphur (see Fig- 
ure |2j. Since this centrally concentrated 56 Ni is enshrouded by a 
much larger mass envelope than the 56 Ni from the He-detonation, 
the outwards diffusion time will be longer. Thus the light curve 
will evolve on longer timescales than for the p-Ia events predicted 



by iBildsten et al.1 d2007t) . IShen & Bildstenl J2009h and IShen et al.l 
d2010l) ; see SectionlBI 



Qualitatively similar results are found for our CSDD-L model. 
In this case, 48 Cr is the dominant radioactive product from burning 
of the He layer. Also, the low density of the CO core means that the 
56 Ni mass produced in the core detonation is now smaller than the 
mass of radioactive elements produced in the He-shell. Neverthe- 
less, the core 56 Ni mass is not negligible and affects the synthetic 
observables, as will be discussed in Section[5] 
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Table 3. Masses of important elements and isotopes in the ejecta for our converging-shock double-detonation (CSDD), edge-lit double-detonation (ELDD) 
and He-only detonation (HeD) simulations for our standard (S) and low-mass (L) models. The upper part of the table lists the masses of material originating 
in the outer He layer while the lower part refers to the products of the CO core. Only the unbound core material is included for the HeD models, having a total 
mass of 0.041 and 0.029 Mq for models HeD-S and HeD-L, respectively. 







CSDD-S 




ELDD-S 




HeD-S 




CSDD-L 




ELDD-L 




HeD-L 




He-layer ejecta 


56 Ni(M ) 


6.5 x 10" 


2 


7.2 x 10- 


2 


6.5 x 10- 


2 


2.8 x 10- 


3 


7.6 x 10- 


-3 


2.8 x 10- 


-3 


composition 


52 Fe (M Q ) 


1.4 x 10" 


2 


1.2 x 10- 


2 


1.4 x 10- 


2 


9.4 x irr 


3 


1.4 x 10" 


2 


9.4 x 10" 


-3 




48 Cr(M ) 


8.6 x 10" 


-3 


7.7 x 10- 


3 


8.6 x 10" 


-3 


1.8 x 10- 


2 


1.9 x 10- 


2 


1.8 x 10- 


2 




Ti (Mq) 


3.7 x 10- 


-3 


3.3 x 10- 


3 


3.7 x 10- 


-3 


1.1 x 10- 


2 


9.6 x 10- 


-3 


1.1 x 10- 


2 




Ca(M ) 


8.8 x 10- 


-3 


8.2 x 10- 


3 


8.8 x 10- 


-3 


2.7 x 10- 


2 


2.5 x 10- 


2 


2.7 x 10- 


2 




S(M ) 


1.7 x 10- 


-3 


2.5 x 10- 


3 


1.7 x 10- 


-3 


6.0 x 10- 


3 


6.5 x 10" 


-3 


6.0 x 10- 


-3 
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1.5 x 10- 


-3 


6.0 x 10- 


3 


1.5 x 10" 


-3 


1.9 x 10- 


3 


4.9 x 10" 


-3 


1.8 x 10" 


-3 




Ma (M,-,) 
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4 
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-4 


5.9 x 10- 


4 


1.3 x 10- 


-3 


5.9 x 10- 
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-3 


1.4 x 10- 


3 
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-3 


6.9 x 10- 


4 


1.4 x 10- 


-3 


6.9 x 10- 


-4 
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1.9 x 10" 


4 


1.5 x 10- 


4 


1.8 x 10- 


-4 


1.5 x 10- 


3 


9.4 x 10" 


4 


1.5 x 10" 


-3 




He (M ) 


8.0 x 10- 


-2 


7.3 x 10- 


2 


8.0 x 10" 


2 


1.2 x 10" 


1 


1.1 x 10- 


1 


1.2 x 10" 


-1 


CO core ejecta 


56 Ni(M Q ) 


1.5 x 10" 


1 
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:s 






1.9 x 10- 


2 


1.5 x 10" 
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composition 


52 Fe(M ) 
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3 
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4.5 x 10" 
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2 
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2 


1.4 x 10" 


2 




He (M ) 


3.5 x 10- 


-4 










1.9 x 10- 


5 











4 ADDITIONAL MODELS 

For comparison of the results obtained with our CSDD models, we 
have performed additional simulations to quantify the observable 
properties of alternative explosion mechanisms. These were set up 
and carried out in an identical manner except that the detonation of 
the CO core was handled differently. 



4.1 He detonation only; HeD models 

Although the simulations described in Section [3] suggest that det- 
onation of the surrounding He layer will trigger a secondary deto- 
nation of the core, the difficulty in det ermining whether a deto na- 
tion is initiated must be recognized (e.g. lSeitenzahl et al.l2009bl see 
also Section RTTV Therefore, as an experiment, we also performed 
calculations in which it is assumed that no core detonation is ignited 
(hereafter, our HeD models, which were carried out for both our S 
and L initial systems). These models are an important compa rison 
point since they are the realisations of the lBildsten et al J ( 120071) p-Ia 
explosion scenario that correspond to our CSDD models. 

For these simulations, the He detonation was ignited exactly 
as before and produces very similar nucleosynthetic yields to the 
He layer in the CSDD model (Table[3}- Since it is assumed that no 
core detonation takes place, most of the underlying CO core is unaf- 
fected by the He detonation and remains tightly bouna[fl This bound 
material is not included in our radiative transfer simulations, which 
involve only the homologously expanding ejecta. However, part of 
the CO core (~0.041 M for HeD-S and 0.029 M for HeD-L) is 



6 We note, however, that even in this case the heating in the shock conver- 
gence leads to burning a small fraction of the mass to 56 Ni (~10 — 4 M 
for our S-model) and intermediate-mass elements (~10~ 3 M ). 



unbound as a result of kinetic energy transf erred from the He det - 
onation to material of the CO core (see also lWooslev et al]|l98ot) . 
Therefore, some material from the core is still present in the ejecta 
and dominates the composition at low velocities (see lower right 
panel of Figure |2). 



4.2 Edge-lit core detonation; ELDD 

An alternative to the CSDD model is that the He detonation di- 
rectly ignites an inward propagating detonation at the edge of the 
CO core. Whether such a detonation can be ignited depends on 
many factors including the density at the edge of the CO core 
(pb), the composition and the geometry. ID simulations have sug- 
gested that it is most likely to happen if the ig nition point is 
some way above the base of the He-shell (see e .g. iNomotoll 19821 ; 
lLivndll99d : lBenzlll997l : rGarcia-Senz et al.lll999h . Edge-lit detona- 
tion has been found in some multi-dimensional simulations (e.g. 
iLivne & Glasnedfl99lt) although i t may be harder when the ID 
symmetry is broken "( Forcadall2007h . 

Edge-lit CO detonation requires that densities of at least 



10 g cm and critical tempera tures of several bil lion Kelvin are 
reached in the outer CO material ( tRopke et alj2007l) . In the simula- 
tions for our S-model, we do find that some regions at the very edge 
of the CO core are heated to temperatures > 2 x 10 9 K and that 
densities in excess of 10 6 g cm -3 are reached in some places (par- 
ticularly close to the He detonation convergence point on the — z- 
axis). However, these hot/dense conditions appear in only a handful 
of our nucleosynthesis tracer particles and it is therefore marginal 
whether critical volumes for detonation are really reached. More- 
over, since our models are not based on evolutionary calculations, 
we cannot predict at what height in the He-layer ignition of the det- 
onation is most likely to occur - we have simply chosen to ignite 
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our He detonations at the base of our He layers. Thus, our sim- 
ulations are ill-suited to determine whether edge-lit detonation is 
probable in our particular systems. 

Nevertheless, we can investigate the observable consequences 
if such an edge-lit detonation were to occur - to do this we per- 
formed simulations in which a CO detonation was ignited by hand 
at the edge of the core immediately below the ignition point of the 
He detonation. The two detonations were ignited simultaneously 
but modelled independently. We will refer to the results of these 
simulation as our edge-lit double-detonation models (ELDD-S and 
ELDD-L, for our two initial systems, respectively). The important 
difference from the CSDD models is that there is no strong shock 
convergence in the core prior to ignition of the CO detonation. This 
means that the densities in the core remain close to their initial val- 
ues, which are too low to lead to significant 56 Ni production in 
the core. Thus, although sufficient energy is released to unbind the 
core, the core is primarily burned to intermediate mass elements 
(see Tableland Figure [2} . In addition, these models have slightly 
more complete burning of the He shell material, a consequence of 
additional heating of the burning region just behind the He detona- 
tion by oblique shocks generated from the CO detonation. 



5 SYNTHETIC OBSERVABLES 
5.1 Light curve morphology 

Figures [3] and [4] show synthetic bolometric (ultraviolet-optical- 
infrared; hereafter UVOIR), optical- and infrared-band light curves 
computed for the three explosion mechanisms (CSDD, ELDD, 
HeD) for our two model systems (S and L, respectively). Here, we 
show the angle-average synthetic light curves - the dependence of 
the light curve properties on observer inclination will be discussed 
briefly in Section |531 

It is apparent from Figure[3]that the CSDD scenario leads to a 
very different transient from a p-Ia model (i.e. the HeD model) for 
ou r S -model. I n agre ement with the calculations for p-Ia models 
by IShen et all fcOld) . our HeD-S bolometric light curve reaches 
peak around 8 days after explosion and then decays fairly rapidly, 
dropping by ~1 mag during a two- week period after maximum. In 
contrast, the CSDD-S model takes several days longer to reach peak 
and remains bright for an extended period - significant bolometric 
decline does not commence until ~30 days after explosion. The 
luminosity of the CSDD-S model is always higher than the HeD- 
S model. Similar conclusions are drawn from our low-mass model 
(compare CSDD-L and HeD-L in Figure [4]l. Here the scale of the 
effect is less extreme but the slower light curve evolution is still 
very apparent in the optical bands. 

These differences can be understood as consequences of 56 Ni- 
rich material produced in the core detonation (see Table [3} . Decay 
of 56 Ni in the core produces a comparable amount of energy to that 
supplied by decay of radioactive nuclei in the outer layer of He- 
burning products. Moreover, the 56 Ni in the core is deep inside the 
ejecta meaning that this energy takes longer to diffuse outward and 
7-ray trapping is more effective for a longer period of time. This 
causes the slow light curve evolution. The scale of these effects is 
large and easily observable, corresponding to differences in excess 
of a magnitude in most bands at post-maximum epochs for our S- 
model. 

The influence of the core material in the ELDD models is 
considerably more subtle. From about 8 days after explosion, the 
ELDD-S bolometric and optical light curves are much more sim- 
ilar to the HeD-S than CSDD-S light curves. This is because the 



ELDD-S light curve is predominantly powered by the radioactive 
nuclei produced in the He detonation; the small mass of 56 Ni in 
the core of the ELDD-S model only becomes a noticeable energy 
source well after maximum light (at times greater than ~20 days 
after explosion, the ELDD-S optical band light curves are system- 
atically brighter than those of the HeD-S model). Compared to the 
HeD-S model, the ELDD-S light curves reach peak slightly ear- 
lier (and are a few tenths of a magnitude brighter at peak). This 
follows from the different velocity distribution of the products of 
He burning in these models. In the ELDD-S model, the He-layer 
ashes have all been pushed out to relatively high velocity by the 
underlying core material (see Figure^. Consequently, the outward 
column density from the surface layer of 56 Ni and 52 Fe is smaller 
in the ELDD-S model, leading to earlier 7-ray escape and down- 
turn of the UVOIR light curves. This effect, however, is relatively 
modest in scale: overall the optical band ELDD-S light curves are 
not very different from those of the HeD-S model. Similar conclu- 
sions can be drawn by comparing the HeD-L and ELDD-L optical 
band light curves in Figure [4] 

An important difference between ELDD and HeD light curves 
manifests in the near-infrared (e.g. the J-band for our S-model 
shown in Figure [3]l- In the near-infrared, a significant fraction of 
the emission in both the CSDD and ELDD models is provided 
by the intermediate-mass elements in the relatively cool and dense 
ejecta from the CO core, particularly around maximum light for the 
CSDD-S and ELDD-S simulations. This emission from the core is 
powered by a combination of energy injected by the radioactive ma- 
terial in the core ejecta (significant for the CSDD-S model) and irra- 
diation by the overlying He-shell ejecta - around maximum light, 
re-radiation of energy originating from the He-detonation ash is 
the dominant source of NIR emission for the ELDD-S simulation. 
Since only a very small mass of material was unbound from the 
CO core in our HeD models, these processes are largely absent in 
the HeD models, making them fainter at these wavelengthfl Thus, 
near-infrared data could be particularly valuable when hunting for 
direct observable signatures of the core detonation. 

5.2 Colours and spectra 

Maximum light colours only weakly discriminate between our 
models. In all cases, B — V, V — R and V — I are positive around 
maximum light (see Figure|5]for our S-model results) and differ by 
at most a few tenths of a magnitude between the explosion scenar- 
ios. 

After maximum light, the B — V colour rapidly becomes more 
positive because of the decline in B-band for all models. Evolution 
of the redder optical colours is complex but qualitatively similar to 
th e colour evolut i on fou nd in the double-detonation models studied 
bv lKromer et al. In particular, our CSDD-S model displays 

similar colour evolu tion to the lowest mass models in figure 3 of 
iKromer et al] d2010h . The ELDD-S and HeD-S models also show 
evolution to redder colours immediately after maximum light. This 
is both faster and more pronounced than in the CSDD-S model and 
both models show extremely red V — I colours within two weeks 



Although the remnant CO WD will remain in the centre of the ejecta, it 
is expected to be too small to intercept a significant fraction of the radiation 
created in the rapidly expanding He detonation ash. Therefore, in contrast 
to the core ejecta in the CSDD and ELDD models, the WD is not expected 
to provide an effective target for re-radiating a significant fraction of the 
emission around maximum light. 
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Figure 3. Angle-averaged light curves for our explosion simulations based on our S-model (CSDD-S, ELDD-S and HeD-S) in bolometric (ultraviolet-optical- 
infrared, UVOIR), B-, V-, R-, I- and J-bands. 
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Figure 4. As Figure|3]but showing results for simulations with our L-model (CSDD-L, ELDD-L and HeD-L). 



of maximum light, a consequence of strong cooling emission by 
the Ca II infrared triplet contributing to the T-bancffl. 



As in the iKromer et alj bold) calculations, the red colours 
in our models are a consequence of effective line-blocking at blue 
wavelengths by iron group elements in the outer ejecta. In particu- 



8 Some caution must be applied to the interpretation of our prediction of 
very powerful Ca II emission at ~30 days after explosion in the ELDD 
and HeD models - by these epochs the ejecta are sufficiently dilute that 
forbidden line emission may contribute significantly to the line cool- 



ing. Such emission is neglected in the current implementation of ARTIS 
iKromer & Sirnll2009t) meaning that the calculations may overestimate the 
strength of Ca II emission at late epochs. 
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Figure 5. Optical colour evolution (B — V , V — R, V — I) for our three explosion simulations (CSDD-S, ELDD-S, HeD-S). 



lar, the He detonation in all models yielded significant masses of Ti 
and Cr that, along with Ca, strongly influence the spectrum. This is 
illustrated in Figure [6] which shows the spectra of the HeD-S and 
HeD-L models at 8 days after explosion (the result is very simi- 
lar for our other models). The colour coding in the figure indicates 
which elements were responsible for the last physical interactions 
of escaping Monte Carlo quanta in our radiative transfer simula- 
tions, making clear the dominance of elements with Z = 20 to 24 
(Ca to Cr) in shaping the emergent spectrum. 

Figure [7] compares the spectra for our explosion scenarios at 
two epochs, 8 and 20 days after explosion. At the earlier epoch 
(around peak brightness for the ELDD and HeD models), the spec- 
tra are all quite similar with strong absorption in the Ca II in- 
frared triplet and Ti II lines (e.g. the characteristic trough around 
4200 A). There are some subtle differences, however - for exam- 
ple, the ELDD calculations (for both S- and L-models) typically 
show higher velocity line features than the HeD simulations, a con- 
sequence of the He-layer ejecta having been pushed to higher ve- 
locity by the CO detonation in the ELDD model. In general, our 
spectra are qualitatively similar to the maxi mum light spectra for 
the low-mass models of lKromer et al.l d2010h but with a noticeable 
reduction in the role played by silicon and sulphur. 

By 20 days, the differences between the models are more ap- 
parent in the spectra. In particular, the CSDD-S model begins to 
show additional, relatively narrow line features (e.g. between 5600 
and 5900 A) -these are formed in the slowly expanding ejecta from 
the CO core detonation. In contrast, the ELDD-S and HeD-S spec- 
tra remain dominated by broad line features that form in the outer 
ejecta. In addition, the ELDD-S (and ELDD-L) spectra now very 
clearly show higher velocity features than the corresponding HeD 
models. 



5.3 Observer inclination 

As discussed by |Finketal1 ( l2010l) . single-spot ignition of the He 
layer leads to a global asymmetry in the explosion that influences 
the ejecta from both the He layer and the core (see Figure |2j- In 
particular, the He-layer detonation ash tends to sweep around the 
CO material (in the same sense as the laterally propagating He det- 
onation), leading to a more geometrically extended layer of He- 
burning products around the pole opposite to the He ignition point. 
This asymmetry affects th e light curves and sp ectra, particularly at 
blue wavelengths (see e.g. lKromer et aljboioh . Figure[8]illustrates 
this for our CSDD-S, ELDD-S and HeD-S models in the UVOIR, 
B-, V- and J-band light curves. The UVOIR curves are sensitive 
to observer inclination by several tenths of a magnitude while the 



bluer optical bands (B- and y-band) are affected more strongly (up 
to ±0.5 mag variation around t he angle-averaged in B-band). As 
in the simulations discussed by iKromer et alj a2010h . the colours 
are bluest and the light curves decline most rapidly when viewed 
from the side on which the He detonation was ignited (i.e. from 
the ±z-direction). At red wavelengths, observer inclination is less 
important and becomes negligible in the near-infrar ed (see lowest 
panels in Figure[8j, in accordance with the findings of lKromer et al.1 

boioh . 

The influence of observer orientation is more complex in the 
CSDD-S model than the HeD-S and ELDD-S models. In that case, 
maximum light in both B- and V-band occurs significantly later 
(and is brighter in V) when viewed from the — z-direction. The 
broad light curve peaks in these models are sustained by radiation 
diffusing out from the 56 Ni-rich inner parts of the ejecta (the CO 
ash). Since the CO detonation is off-centre (displaced down the — z- 
direction; see Figure|2]and Table|2j, more of this radiation emerges 
in the — z-direction, leading to brig hter extended op tical emission 
when viewed from this direction (cf. lSim et alj|2007h . 

In the L-model simulations, the burning of the He-layer is sig- 
nificantly more complete in the ash around the — z-axis than in 
other directions (the 56 Ni mass fraction is small in most regions 
of the the He-layer ash but it becomes significant around the — z- 
axis because of the enhanced burning in the region where the He 
detonation converge a. Nevertheless, the effect of orientation on 
the L-model light curves is qualitatively similar to that found in the 
S-model and the scale of variation is comparable (up to ±0.5 mag 
in the blue bands and negligible in the near-infrared). 

Overall, the influence of observer inclination is modest com- 
pared to the difference in light curve morphology between the 
CSDD and ELDD/HeD models. In particular, by ~20 days after 
explosion, the light curves are fairly orientation-independent in all 
bands. Thus, observer orientation should not severely hinder obser- 
vational discrimination between the CSDD and ELDD/HeD sce- 
narios. However, it does affect the light curve shapes on a sim- 
ilar scale to the difference between the ELDD and HeD models 
and therefore will complicate attempts to distinguish these mecha- 
nisms. 



9 Due to reduction of surface area, the detonation shock strengthe ns and 
the d etonation becomes over-driven around the convergence point jLivnel 
Il997l) . This affects the yields since the densities are just below the critical 
densities at which 56 Ni is produced in a detonation. 
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Figure 6. Optical spectra of model HeD-S (left) and HeD-L (right) at 8 days after explosion. The upper boundary of the coloured loci along the bottom of the 
plots are the synthetic spectra. The colour coding under the spectra identifies the elements with which escaping Monte Carlo quanta in each wavelength bin 
last interacted in our radiative transfer simulation. The coloured region along the top of the plots indicates which elements were last responsible for removing 
Monte Carlo packets from a particular wavelength bin (specifically, it shows the distribution of photon wavelengths that escaping packets had prior to their 
last interaction). The colour bar on the right indicates the colour coding used for each atomic number (Z). 



6 DISCUSSION AND CONCLUSIONS 



6.1 Core detonation 



In previous work, we have studied the double-detonation model ap- 
plied to sub-Chandrasekhar mass CO white dwarfs with masses 
~1 Mq: this scenario predicts transi ents with a range of bright- 
ness compatible with SNe la (e.g. Hofli ch & Khokhlovj 1 19961 : 
iNugent et all 1 1991 iKromer et al.ll20ld : IWooslev & Kasenl Y20\% 
and appears to be a ro bust explosion mechanism iLivne & Arnetj 
ll995l : lFink et alj201(t) . 

Here we have investigated the plausibility of double- 
detonation models for systems with relatively low CO core mass 
and hi gh He-layer mass, following the methodology of iFink et al.l 
d2010l) . We concluded that, if detonation of the He layer occurs, 
the resulting compression and heating of the CO core by inwardly 
propagating shocks can produce sufficiently high densities and tem- 
peratures that core detonation may occur, even for core masses as 
low as 0.45 Mq . Tog ether with previous results Jlivne & Arnetj 
ll995UFinketaU20ia) this suggests that detonation of a n accreted 
He layer (as in the p-Ia scenario of iBildsten eTai]|2007l) could be 
accompanied by explosion of the underlying core for all CO core 
masses that are commonly realized in nature. 

Although our results suggest that core detonation is probable 
in all cases, they do not prove that it must always occur. For ex- 
ample, strong rotation might inhibit the shock convergence. Alter- 
natively, the converging shocks might heat the material around the 
putative detonation point sufficiently that burning occurs prior to 
ignition of a detonation. Depending on the geometry, this might 
mean that the CO hot spot is completely enshrouded in nuclear ash, 
starving the detonation of fuel such that it is not able to propagate 
and incinerate the whole star. To investigate these possibilities in 
the future will require high-resolution, three-dimensional hydrody- 
namical/nucleosynthesis simulations. It will also be important to 
study the double-detonation mechanism for systems with ONe WD 
cores - in this case it will be harder to ignite a core detonation, 



potentially opening an alternative parameter-space for He-layer ex- 
plosions. 



6.2 Observable signatures of core detonation 

For large CO core masses (>0.9 Mo), the double-detonation 
model predicts light curves that are primarily powered by the large 
mass of 56 Ni synthesised in the core detonation. Therefore these 
events will look very different from predictions of the p-Ia sce- 
nario for explosion of the corresponding He layer without the core 
detonation. For lower mass CO cores, however, the o6 Ni mass pro- 
duced in the core explosion is reduced such that some portion of the 
light curves will be predominantly powered by the radioactive nu- 
clei formed in the burning of the He layer, as in the p-Ia model. To 
study this, we performed radiative transfer simulations for three ex- 
plosion scenarios: He detonation followed by core detonation trig- 
gered following shock convergence (CSDD models), He detona- 
tion followed by prompt edge-lit core detonation (ELDD models) 
and He detonation with no core detonation (HeD models). We in- 
vestigated these explosion scenarios for two initial model systems 
(M C o = 0.58 and 0.45 M together with He layers of 0.21 M ). 

We found that the CSDD models can be easily distinguished 
from the ELDD and HeD models by their light curve morphologies: 
compression of the core prior to detonation means that the CSDD 
model yields relatively large masses of 56 Ni from the core (for both 
our initial systems, the mass of 56 Ni produced in the core is com- 
parable to the total mass of radioactive nuclei produced in burning 
of the He layer). This central 56 Ni causes the light curves to fade 
more slowly than in the other models. Consequently, this scenario 
can be easily distinguished observationally from an explosion of a 
He surface layer alone. 

In contrast, our ELDD models show that if core detonation 
is ignited without significant pre-compression, its influence on the 
observables is much more subtle for the low-Mco systems we con- 
sider. The modest differences between the synthetic observables for 
our ELDD and HeD models can mostly be attributed to the differ- 
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Figure 7. Optical spectra for our S- and L-models (top and bottom, respectively; CSDD in black, ELDD in red and HeD in blue) at 8 and 20 days after 
explosion (left and right panels, respectively). The flux scale is normalised to the peak of the CSDD model in all panels. 



ing structure of the ejecta layers produced in the He detonation. In 
particular, explosion of the core in the ELDD model clears out the 
He detonation products from the inner ejecta leading to light curves 
that evolve a little more rapidly and spectroscopic features that are 
broader and more blue-shifted. Thus rapidly evolving thermonu- 
clear transients can be produced by edge-lit double-detonation of 
low-mass systems. Distinguishing them from pure He-layer deto- 
nations could best be done via the evolution of spectral line shapes 
and infrared photometry (see Section[5}. 

It is noteworthy that the difference between our CSDD and 
ELDD models is stronger for our more massive initial system 
(Mtot = 0.79 Mq) than our extremely low-mass system (M tot = 
0.66 Mq). In both cases, the pre-explosion CO densities are too low 
to yield large 56 Ni-masses in a prompt detonation. However, for 
our more massive system, the compression by converging shocks 
is able to raise the density sufficiently to produce Ni in a fairly 
substantial fraction of the core (our CSDD-S model) - this leads 
to the differences in shape between our CSDD-S and ELDD-S 
light curves. For the low-mass system, the same effect occurs but 
is less dramatic because the initial CO densities are so low that 



even with shock compression only a small fraction of the core is 
burned to 56 Ni. For even more massive systems (M to t ^ 0.9 Mq), a 
large fraction of the CO fuel will be at sufficiently high densities to 
produce 56 Ni re gardless of shock compression (see e.g. table 1 of 
ISim et alj2dToh . Double detonations of such systems will therefore 
always produce significant masses of 56 Ni from the core, meaning 
that the differences between the ELDD and CSDD scenarios should 
be relatively small for massive systems. Thus the best opportunity 
to observationally distinguish between the ELDD and CSDD mech- 
anisms will be in systems where the pre-explosion densities in the 
core are close to but below the critical densities at which o6 Ni is 
produced in a detonation. 



6.3 Relation to known transients 

The goal of this study has been to investigate whether secondary 
core detonation is likely for low mass CO cores and to predict the 
influence of CO core detonation on synthetic observables. We have 
not yet conducted an exploration of parameter space as required 
to attempt a quantitative comparison with observations. We can, 
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Figure 8. Light curves for selected observer inclinations to our CSDD-S, ELDD-S and HeD-S explosion models (left to right, respectively) in UVOIR, B-, 
V- and J-bands. In each panel three light curves are shown for different observer orientations: viewed down the equator (black), from the +,z-direction (blue) 
and the — ^-direction (red). In all models the He ignition spot was on the 4-z-axis and the angle-averaged light curves (see Figure 3) are very similar to those 
for an equatorial line of sight. 
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however, comment qualitatively on the relation of our synthetic ob- 
servables to the properties of known astrophysical transients. For 
reference, we tabulate important light curve shape parameters (rise 
times and decline timescales) for our simulations in Table|4] 

IPerets et alj d2010h proposed that the Type lb event SN 2005E 
could be attributed to an explosion of an accreted He-layer on a 
WD. Its peak brightness suggests that only a few thousandths of a 
solar mass of 56 Ni, 52 Fe and 48 Cr were produced, if such a model 
is applicable. A p-Ia like exp losion was also di scussed as one possi- 
bility to explain SN 2008ha l lFolev et alj2009h . this event was even 
fainter than SN 2005E and spectroscopically different, showing low 
expansion velocities and no evidence of He. In both cases, however, 
the rise time was estimated to be < 10 days and the light curve 
decline parameter was AM^ ~ 2 mag. These rapid timescales 
clearl y invite proper c o mpari son with models for He-layer explo- 
sions. IWaldman et all ( feOllri attempted to model SN 2005E in 
the context of the p-Ia scenario. They investigated a model with 
Mco = 0.45 M and M He = 0.2 M (their C0.45HE.2 model) 
and found that such a model might be able to account for the 
peak brightness of SN 2005E. However, they also found that, al- 
though significantly faster than normal SNe la, the decline rate of 
SN 2005E was slower than could be easily explained with a p-Ia 
model. 

SN 2005E l lPerets et all |201Ch was significantly fainter 
(MboI > — 15 mag at peak) than our model s. Although the param- 
eters of our L-system are very similar to the lWaldman et al.1 fcOllh 
C0.45HE.2 simulation, our explosions are significantly brighter. 
The large difference in brightness stems from the fact that 44 Ti is 
the most abundant radioactive nucleus produced in the C0.45HE.2 
simulation while our HeD-L model predicts 48 Cr to be dominant. 
Although only a modest shift in the mean nucleosynthetic yields, 
this strongly affects the brightness since the half-life of 44 Ti is 
much longer than that of 48 Cr. Modification of our nucleosynthesis 
to produce less complete burning (i.e. a lower mass of 48 Cr) could 
be achieved by re ducing the density of the He layer, although this i s 
likely unphysical dShen & Bildstenl 20091 : Iwooslev & Kasenl201ll) . 
Alternatively, significantly polluting the He layer wit h a heavier el- 
emen t such as carbon will a lter the nucleosynthesis l lKromer et al.1 
120101 : IWaldman et al.|[20lTh . For a sufficiently large initial mass 
fraction, adding carbon reduces the typical atomic weight of the 
burning products and leads to more intermediate-mass elements 
rather than iron-group material. Reduced yields of radioactive iron- 
group elements would make our models fainter and alter the spec- 
tral features. We have not investigated such possibilities here. Nev- 
ertheless, our simulations have relevance to the study of SN 2005E 
and similar faint and fast transients, potentially including e.g. 
SN 2008ha - in particular, they show that if a He layer detonation 
in a low-mass system is followed by a core detonation triggered 
by converging shocks, the core material can significantly retard the 
post-maximum decline of the light curves. For example, in our p- 
Ia-like HeD-L model, the B-band decline during the 15 day period 
after maximum light is AMf 5 = 3.4 mas while for our CSDD-L 
model, AM^ = 2.5 mag (the scale of this effect is even more 
dramatic in our S-model simulations; see Table [4]l . Thus, in future 
studies, it will be important to consider whether CSDD explosions 
of low-mass systems might be able to account for faint thermonu- 
clear transients whose light curves decline more slowly than can be 
explained by p-Ia models. 

It has also been suggested that p-Ia explosions may ac- 
count for a class of very rapi dly evolving thermonu clear explo- 
sions that includes SN 2002bj jPoznanski et"ai]|2010l) SN 2010X 
dKasliwal et all |201Ch . SN 1939B and SN 1885A dPerets et all 



l201lh . Thes e events are consider ably brighter than SN 2005E: 
for example. iKasliwal et all d2010h estimate M r ~ —17 mag for 
SN 2010X while SN 2002bj is around 1.5 mag brighter. Moreover, 
the light curves of these events have a very rapid post- maximum de- 
cline (for example, AM§ « 3.2 mag for SN 2002bi; |p"erets et al.1 
I20T l]). The rapid post-maximum decline of this class of object is in- 
consistent with our CSDD simulations, particularly for our brighter 
S-model. However, it may be compatible with the rapid decline 
in our ELDD (or HeD) models (see Table |4). Therefore focused 
modelling may be warranted to properly investigate whether ELDD 
models could be applicable to this class of transients. 

6.4 Future work 

There remain many open questions to be addressed. Most impor- 
tant, perhaps, is detailed study of burning of the He layer. The He 
burning products play a critical role in determining the observa- 
tional properties of double-detonation explosion models but their 
yields are sensitive to many i ssues including the exact densities and 



composition of the He layer dShen & Bildste n 2009; Kromer et al.1 
l2010l : IWooslev & Kasenll201ll ; IWaldman et aliboilh . th e direction 
in which the detonation propag ates (radial or azimuthak lFink et"al] 



n pri 

l2010l : IWoosiev" & Kasen 20 1 1|) and the stru cture of the detonation 
front. For example. lKromer et al. and lWaldman et al I fcOllh 

both showed how introducing C to the He-layer prior to burning can 
alter the final composition - for double detonation of models with 
massive CO cores, this dramatically affects the colours and can lead 
to much improved agreement with observed SNe la spectra. 

In addition, we have shown that the observational conse- 
quences of core detonation in low-mass CO cores are quite dif- 
ferent if the core detonation is edge-lit rather than triggered by 
shock compression. Although the triggering of the secondary deto- 
nation in our CSDD models is based o n previous studies of the nec- 
essary conditions for CO detonation dNiemever & Wooslevll 19971 : 
iRopke et al.l2007l) , we have not investigated the physical plausibil- 
ity of edge-lit detonation. Further study of this is clearly warranted. 
This is of primary interest for models with low CO core masses 
(e.g . the S-model studie d here rather than those of lFink et alfeoiol 
and iKromer et al. I I2010I) since the difference between the CSDD 
and ELDD mechanisms will be largest in this case (see Section|5]l- 
However, it does have some relevance to the study of more massive 
cores since the position of an off-ce ntre ignition imprints d istinctive 
signatures on the explosion ejecta dChamulak et al.l201 lh . 

It will also be important to quantify the relative frequency of 
bright (i.e. SNe la-like luminosity) and fainter (i.e. p-Ia-like lu- 
minosity) explosions as a check on the plausibility of the double- 
detonation scenario contributing to both populations. The core/He- 
layer mass combinations used in this study are close to the low- 
mass extreme for which potential progenitor systems can be re- 
alised in nature. For low core masses, the minimu m mass of ac- 
creted He required for detonat i on be comes large dBildsten et al.1 
120071) . As noted by IShen et al.1 d2010l) . systems with massive He 
layers ( > 0.1 M©) are not expected to be reached via the evolution 
of systems in which a CO WD accretes from a He WD donor. Our 
models, however, might be realised in sy s tems where the donor i s 
a He-burning star dlben & Tutukovlll99ll ; lshen & Bildstei]|2009b. 
When investigating potential progenitors for SNe la. iRuiter et al.1 
d201 lh found that the population of potential double-detonation 
systems with He-burning star donors is sub-dominant but can pro- 
vide a significant event rate in young stellar populations. If ex- 
tended to less massive CO cores, such population synthesis stud- 
ies could estimate the relative occurrence of progenitors for the 
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Table 4. Light curve rise times and decline rate parameters for the models. tuvoiR, max and tg m . dx are the times after explosion to maximum light in UVOIR 
and B-band, respectively. AM^ 0,R and AMK parametrize the light curve decline rate in UVOIR and B-band light (specifically, they give the increase in 
magnitude during the 15 days after maximum light in the appropriate light curve). 



Parameter 


CSDD-S 


ELDD-S 


HeD-S 


CSDD-L 


ELDD-L 


HeD-L 


tuvoiR, max (days) 


12 


7.1 


8.7 


11 


7.6 


9.9 


AM UVOIR (mag) 


0.4 


1.8 


1.8 


1.3 


2.0 


1.9 


*B,max (days) 


8.1 


6.0 


7.7 


5.4 


5.2 


6.4 


AMg (mag) 


1.2 


3.2 


3.1 


2.5 


4.0 


3.4 



double-detonation scenario leading to both bright and fainter ther- 
monuclear transients, a prediction that can be tested by the current 
and future generations of wide-field transient surveys. 



ACKNOWLEDGMENTS 

This work was supported by the NCI National Facility at the Aus- 
tralian National University, the Deutsche Forschungsgemeinschaft 
via the Transregional Collaborative Research Center TRR 33 "The 
Dark Universe", the Excellence Cluster EXC153 "Origin and 
Structure of the Universe" and the Emmy Noether Program (RO 
3676/1-1). Parts of the simulations were carried out at the John 
von Neumann Institute for Computing (NIC) in Jiilich, Germany 
(project HMU14). 

S AS thanks Ivo Seitenzahl, Ken Shen, Lars Bildsten and Brian 
Schmidt for illuminating conversations. We thank Riidiger Pakmor 
for developing tools to map nucleosynthesis tracer particle yields 
for our explosion models. 

We thank the referee for useful comments that helped clarify 
the presentation of our work. 



REFERENCES 

Benz W., 1997, in Ruiz-Lapuente P., Canal R., Isern J., eds. Ther- 
monuclear Supernovae Vol. 486 of NATO ASIC Proc, Three- 
dimensional simulations of core ignition in sub-Chandrasekhar 
mass models. Kluwer Academic Publishers, Dordrecht, pp 457- 
474 

Bildsten L., Shen K. J., Weinberg N. N., Nelemans G, 2007, 
ApJL, 662, L95 

Chamulak D. A., Meakin C. A., Seitenzahl I. R., Truran J. W., 

2011, ArXiv e-prints 
Dgani R., Livio M., 1990, ApJ, 361, 540 
Fink M., Hillebrandt W., Ropke F. K., 2007, A&A, 476, 1 133 
Fink M., Ropke F. K., Hillebrandt W., Seitenzahl I. R., Sim S. A., 

KromerM., 2010, A&A, 514, A5 3 
Foley R. J. et al., 2009, AJ, 138, 376 

Forcada R., 2007, in Supernovae: lights in the darkness Edge-lit 
double detonation in subchandrasekhar-mass models for type ia 
supernovae 

Forcada R., Garcia-Senz D., Jose J., 2006, in International Sym- 
posium on Nuclear Astrophysics - Nuclei in the Cosmos Single 
point off-center helium ignitions as origin of some Type Ia su- 
pernovae 

Garcfa-Senz D., Bravo E., Woosley S. E., 1999, A&A, 349, 177 
Hillebrandt W., Niemeyer J. C, 2000, Ann. Rev. Astron. Astro- 

phys., 38, 191 
Hoflich P., Khokhlov A., 1996, ApJ, 457, 500 
Iben Jr. I., Tutukov A. V., 1991, ApJ, 370, 615 



Kasliwal M. M. et al., 2010, ApJL, 723, L98 

Keller S. C. et al., 2007, PASA, 24, 1 

Kromer M., Sim S. A., 2009, MNRAS, 398, 1809 

Kromer M., Sim S. A., Fink M., Ropke F. K, Seitenzahl I. R., 

Hillebrandt W., 2010, ApJ, 719, 1067 
Kurucz R. L., Bell B., 1995, Atomic line list 
Livne E., 1990, ApJL, 354, L53 

Livne E., 1997, in Ruiz-Lapuente P., Canal R., Isern J., eds, NATO 
ASIC Proc. 486: Thermonuclear Supernovae Helium and carbon 
detonations in sub-Chandrasekhar white dwarfs, pp 425-440 

Livne E., Arnett D., 1995, ApJ, 452, 62 

Livne E., Glasner A. S., 1990, ApJ, 361, 244 

Livne E„ Glasner A. S„ 1991, ApJ, 370, 272 

LSST Science Collaborations, 2009, ArXiv e-prints: 0912.0201 

Niemeyer J. C, Woosley S. E„ 1997, ApJ, 475, 740 

Nomoto K, 1980, Space Science Reviews, 27, 563 

Nomoto K, 1982, ApJ, 257, 780 

Nugent P., Baron E., Branch D., Fisher A., Hauschildt P. H, 1997, 
ApJ, 485, 812 

Perets H. B., Badenes C, Arcavi I., Simon J. D., Gal-yam A., 

2011, ApJ, 730, 89 
Perets H. B. et al., 2010, Nature, 465, 322 
Poznanski D. et al., 2010, Science, 327, 58 
Rauscher T, Thielemann F.-K., 2000, Atomic Data and Nuclear 

Data Tables, 75, 1 
Ropke F. K, Woosley S. E., Hillebrandt W., 2007, ApJ, 660, 1344 
Ruiter A. J., Belczynski K., Sim S. A., Hillebrandt W., Fryer C. L., 

Fink M., Kromer M., 2011, MNRAS, 417, 408 
Seitenzahl I. R., Meakin C. A., Lamb D. Q., Truran J. W., 2009a, 

ApJ, 700, 642 

Seitenzahl I. R., Meakin C. A., Townsley D. M., Lamb D. Q., 

Truran J. W., 2009b, ApJ, 696, 515 
Seitenzahl I. R., Ropke F. K, FinkM., Pakmor R., 2010, MNRAS, 

407, 2297 

Shen K. J., Bildsten L., 2009, ApJ, 699, 1365 

Shen K. J., Kasen D., Weinberg N. N., Bildsten L., Scannapieco 

E.,2010, ApJ, 715, 767 
Sim S. A., 2007, MNRAS, 375, 154 

Sim S. A., Ropke F. K., Hillebrandt W., Kromer M., Pakmor R., 
Fink M., Ruiter A. J., Seitenzahl I. R., 2010, ApJL, 714, L52 

Sim S. A., Sauer D. N., Ropke F. K, Hillebrandt W., 2007, MN- 
RAS, 378, 2 

Sullivan M. et al„ 201 1, ApJ, 732, 118 

Travaglio C, Hillebrandt W., Reinecke M., Thielemann F.-K., 

2004, A&A, 425, 1029 
Waldman R., Sauer D., Livne E., Perets H., Glasner A., Mazzali 

P., Truran J. W., Gal- Yam A., 201 1, ApJ, 738, 21 
Woosley S. E., Kasen D., 201 1, ApJ, 734, 38 
Woosley S. E., Taam R. E., Weaver T. A., 1986, ApJ, 301, 601 
Woosley S. E., Weaver T. A., 1994, ApJ, 423, 371 



15 

Woosley S. E., Weaver T. A., Taam R. E., 1980, in J. C. Wheeler 
ed., Texas Workshop on Type I Supernovae Models for Type I 
supernovae. pp 96-1 12 



